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Abstract. The classical delta filters used in the current resonant bar experiments 
for detecting GW bursts are viable when the bandwidth of resonant bars is few Hz. 
In that case, the incoming GW burst is likely to be viewed as an impulsive signal in 
a very narrow frequency window. After making improvements in the read-out with 
new transducers and high sensitivity dc-SQUID, the Explorer-Nautilus have improved 
the bandwidth (~ 20 Hz) at the sensitivity level of W 20 /y/Wz. Thus, it is necessary 
to reassess this assumption of delta-like signals while building filters in the resonant 
bars as the filtered output crucially depends on the shape of the waveform. This is 
presented with an example of GW signals - stellar quasi-normal modes, by estimating 
the loss in SNR and the error in the timing, when the GW signal is filtered with the 
delta filter as compared to the optimal filter. 
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1. Introduction 

Till date, the burst data analysis in the narrow-band resonant detectors like Explorer- 
Nautilus is based on the assumption that short gravitational wave(GW) bursts (duration 
of few millisecond to fraction of a second) appear as delta-like signals in the detector 
bandwidth (BW). Such an assumption was made mainly because the short GW bursts 
emit waveforms of unknown shape and further due to detector's narrow band, one could 
safely assume that the signal emits a flat spectrum in the detector BW. However, recent 
improvements in resonant bar detectors, mainly the read-out with new transducers and 
high coupling dc-SQUID pQ, have improved the BW (~ 20 Hz) at the sensitivity level of 
1CT 20 / V Hz, see Fig[TJa). Thus, it is important to reassess the above assumption which 
we demonstrate in this paper with an example of stellar quasi-normal modes (QNM). 

Astrophysical inputs indicate [2j that various physical processes can excite stellar 
QNM during its evolutionary phases, emitting GW in the BW of the resonant bars. 
Such GW may last for a fraction of a second to few seconds. For example, after the SN 
core collapse, during the cooling phase, the proto-NS emits GW as a damped sinusoid 
signal evolving in frequency as well as damping time and can also chirp in the resonant 
bars PHI- In this case, if the data is filtered through a filter matched to a delta-like 
signal - a delta filter - rather than a proper matched filter, the error in arrival timing 
as well as loss in SNR can arise. Here, as a preliminary study, we discuss these issues 
for a simple case of stellar QNM emitting GW h(t) modeled as a damped sinusoid with 
a fixed frequency fo = ujq/2h and damping time r, as given below 



Here, 8(t — t ) implies h(t) ^ for t > t Q (t is the time of arrival) and zero otherwise. 
The damping time r depends on the underlying physical process. 

To carry out the delta filter assessment, we assume the signal parameters f and r 
are known and compare the outputs of the delta filter to that of the matched filter. 

2. System response 

The incoming h(t) can excite the first longitudinal mode of the resonant bar. A low mass 
electrical transducer is attached to the bar to convert this displacement of the bar end 
face into an electrical signal which is further amplified with the SQUID amplifier. Such 
a bar-transducer system is modeled as a coupled harmonic oscillator with two resonant 
modes [3]. The h(t) provides an external force f x (t) = m x Lh(t)/2 to the bar [m x : the 
reduced mass of the system = Mb ar /2, L: the effective length = AL^/n 2 }. 

The electrical output of the transducer is proportional to the relative displacements 
of the transducer and the bar, u(t) = y(t) — x(t) from their equilibrium positions. In 
Fourier domain 0], the response u(t) to the external force f x (t) is obtained by 



h{t) = h sin[cj (t - to)] e~ ( '-' o)/r 0{t - t ) . 



(1) 



U{juj) = W ux (ju)F x {juj) , 



(2) 



| Possible candidates for coincidences with interferometers like VIRGO. 
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where W ux (ju>) is the system transfer function from the input force to the output 
displacement as given defined in Eq.(1.8) of [3]. The F x (ju) [the FT of f x (t)] for h(t) 
defined in Eq. ((TJ) is 

F M = —j—B{ju,) = (1+ ^ r)2 + r2 ^ ex P(-^o) (3) 

where H{juo) is the FT of h(t). The \F x (ju)\ is a Lorentzian which becomes narrower 
(broader), and in(de)creases in height as r in(de)creases, see Fig. [Ha). Conveniently, 
we choose to = 0. The relative displacement of the transducer u(t) is in units of length. 
The electrical signal at the output of the SQUID amplifier is then given by 

v(t) = Ku{t) = M v u (t) in units of V , (4) 

where M v = KM U and M u is the maximum value of u(t) §. 

The expected power spectral density (PSD) of the noise n{t) at the output of the 
electrical chain is given by 0j 

S t (u) = S n + ot 2 S fx \W ux {ju)\ 2 + a 2 S fy \W uy (juj)\ 2 . (5) 

The S n is the broad-band noise contribution from the SQUID and the electrical chain 
with a flat spectrum in the BW and the rest is the narrow-band noise contribution - 
due to the thermal noise of the two mechanical oscillators. The Sf x and Sf y is the total 
noise force spectra due to the Nyquist and the back-action force. The W ux and W uy are 
the system transfer functions as defined in jlj ||. 

3. Matched filter 

Signal detection problem involves computing a statistic - a functional of the data z(t) - 
which when passed through the threshold allows to make the decision of either presence 
[i.e. z(t) = n(t) + v(t)} or absence of signal [i.e. z(t) = n(t)} in the data. 

For known signal in Gaussian-stationary noise, matched filtering is the optimal 
filtering. The matched filter output is given by 

1 r°° 

o(t)=<z,q>=— Z(juj)Q{juj)exp(jujt)duj . (6) 

Z7T J-oo 

The matched filter transfer function for v(t) is Q{jui) = NJJq^u) / 'S t (uj) with Uq{juj) 
as the FT of uo(t). The normalization N u = \Uo(juj)\ 2 / St^du;} -1 is such that 

Max < Mo, q >= 1. In no noise case, z(t) = v(t) and 

o(t)=<v,q>=— — exp(jut)du. (7) 

Z7T J-oo O t [UJ) 

Thus, using the above normalization, we obtain Max[o(t)] = M v , from which one can 
estimate the strength of the input GW ho- The matched filter SNR can be evaluated 
as SNR 2 = (Max[o(t)]) 2 /Var(< n,q >). The variance of the filtered noise is given by 

Var(< n,q>) = E{< n, q > 2 } - (E{< n, q >}) 2 = N u . (8) 

§ K is equal to aAB of Ref. 0] which includes SQUID amplification and transducer constants. 
|| Besides, the other spurious unknown noise sources are treated while filtering the data. 
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Thus, the matched filter SNR is given by 

SNRJ, = M = 1 ffiMi,, (9) 



N u 2ir J-oo S t (uj) 

where V(ju) is the FT of v(t). In terms of the signal and system poles Pi, i — 1, . . . 6 

pi = -w+ + j/r+, p 2 = -p*, P3 = -u-+j/r-, pa = -pi, p 5 = -uj + j/T, p 6 = -p* 5 , (10) 

and the decay times t+,t_ pertaining to the response of the two modes [i.e. f±, 
u± = 27r/±] at the output of the delta filter + , the SNRm is 

h 2 L 2 co 2 K 2 , 
16&, 



SN R^ = "° 



r +Pi , T -Pl , r Ps 



w+ ni,^i, 2 (Pi - ^- rii,^3,4(P3 - pi) rii,^5,6(pi - pD 



4. Explorer-Nautilus Delta filter 

As stated earlier, the delta-filtering is the most natural approach for detecting unknown 
short GW bursts in the narrow-band detector. The delta filters are developed as follows. 

The normalized system response to a delta-like signal (used to build the delta- 
filter) is u 5 (t) such that U 5 {ju) = W ux /M 5 [4 and Max[u s (t)] = 1. Then, the 
delta filter, constructed from this response has the transfer function Q 5 (jw) = 
N s U s *(juj)/S t (u)) with the filter normalization N s = JZo \U S {juj)\ 2 / S t {uj)duj}- 1 . This 
filter construction is such that if an impulse is incident on the bar with v{t) = M v u s (t) 
then the maximum of the filtered output Max(< v, q 5 >) is M v . 

However, when the response of the detector to the damped sinusoid is filtered 
through the delta filter, the filtered output becomes 

o(t) = <v,q >=— — — exp(jwt)dw, 

Z7T J-oo <~>i(k> ) 

—h LuJoKNs f°° cu 4 exp(jut)du 



AnS n M s J-oo (to 2 - pi) (uj 2 - pi) (uj 2 - pi) (uj 2 - p\) (uj - p 5 ) (u - p 6 ) 
We solve this integration by applying the residue theorem and obtain 

. -h Luj KN s T + e-*^+ p\e~^ 
°( f ) = o~o~~h7 ^t" 



8S n M S UJ + ni=3,4(Pl-^)Ilfc=5,6(Pl-PA) 

T _ e -t/r- pie-^- 1 je- t/T v \ e -3^t 
— — V- — (11) 

W- rii=l,2(P3 - Pi) Uk=Sfi(P3 ~ Pk) ^0 rii=l,..,4(p| - Pi) 

The variance of the filtered noise is Var(< n, q 5 >) = N$. Thus, 

M 2 

SNR2 = -^(Max< Mo ,g 5 >) 2 . (12) 
5. Comparison and Numerical plots 

To illustrate, we use the parameters pertaining to Explorer [Feb 2005]; two resonant 
frequencies /_ = 904.7 Hz, /+ = 927.452 Hz, r+ ~ 140 ms, r_ ~ 150 ms, K ~ 
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Figure 1. (a) vs / for r = 50, 250, 450, 650 ms, the higher peak corresponds 

to higher r. (b) Explorer : Two sided PSD in per VSz, (c) SNR m vs t for h ~ 10~ 20 . 
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1.66 x 10 13 V/m and S n ~ 10~ 8 V 2 /Hz. FigHfb) shows the noise PSD of Explorer. The 
BW of the detector at the level of \fS^ ~ 10~ 20 / V Hz is ~ 20 Hz and we define the 
sensitive frequency band as FB : {900, 932}. We note that at / ~ 915 Hz, the sensitivity 
is the worst within FB. We divide our study in 2 cases, (a)/o G FB, (b)/o far from FB. 

FigJTfc) shows the plot of SNRm vs t for fixed ho ~ 10~ 20 . As r increases SNRm 
increases as the signal spends more time in the detector. However, for fo close to f±, 
this increase is sharp as the incoming h(t) excites the resonances and gives more and 
more energy to f± as r increases, see Fig^a). For a given r, the difference in SNR's 
of two incoming GW with frequencies /_ and 915 Hz is related to the corresponding 
S h (f-) and S h (Q15) [see FigIH»]. In FigH^c), the SNR for f = 900 Hz and /„ = 915 Hz 
are similar as the detector sensitivity is similar at those frequencies. For f Q away from 
the resonance, the detector band falls in the tail of the signal Lorentzian giving small 
power to the resonances even at high r. Thus, the increase in SNR is very slow in such 
case. The similar plot can be obtained fixing the signal energy instead of ho- However, 
the signal energy of QNM is itself a function of r. Thus, for clear demonstration of 
dependence of SNR on r, we fix h here. 

5.1. SNR A/ vs SNR 5 

To validate the delta filter, we study the loss in SNR when the signal is filtered with 
the delta filter, i.e. the ratio between Eq. (fT2~]) and Eq.(j0J). In Fig. we plot this ratio 
for case (a) and (b) respectively. We show that when / G FB, the SNR^ is comparable 
to that of the SNRm (assuming the SNR loss of ~ 15%) for all r < 50 ms. This loss 
increases as r increases. Thus, when fo G FB, for high r, the delta filter is far from 
optimal. Contrarily, when f is far from FB, the tail of the signal Lorentzian gives 
relatively flat spectrum in detector BW (for small values of r), similar to a delta-like 
signal. Consequently, delta filter matches the signal giving SNR M comparable to that 

% For detailed discussion, see [1] 

+ In 4 , t + ,t_ are indicated by T3 + ,T3_ respectively, 
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of SNR5 for t as large as ~ 200 ms. However, as r increases, due to the nature of signal 
Lorentzian, the energy given to both the resonances is not the same. This results in the 
decrease in SNR ratio as r increases even for case (b). This loss is related to the error 
in the estimation of h . However, we note that in case(b), the SNR M is well below the 
SNRm obtained for case(a) for the same ho- 

Figure 2. SNR loss in (a) Case (a), (b) Case (b), (c) Error in timing 



(a) £ : [900,915] Hz (b) f„ far from FB (c) 
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5.2. Arrival timing error 

In Fig. El we plot the output of the delta as well as matched filter for case (a) with 
fo = /_ and r = 0.01, 0.2, 1 s. We show in Fig. |3] A (a), when r < t±, the decay of the 
matched filtered output is dominated by the r±(~ 140ms). However, as r increases and 
is > t± [see Fig El A (b) and (c)], the decay time of the filtered output is dominated by 
the signal term which can increase the arrival timing error. However, it is worth noting 
that with the noise, this error not only depends on the decay time of the filtered output 
but also on the SNR. As r increases, the decay time increases but at the same time the 
SNRm also increases (which we wish to investigate in future with simulations). 

The output of the delta filter is shown in Fig|3]B. The output is asymmetric about 
t = to = as the delta filter is causal [as it is not properly matched] unlike the matched 
filter. Mathematically, it can be seen by the the relative sign difference between the 
signal term and the two resonances in Eq. (|TT| . In this case, the filtered output becomes 
maximum when t ~ to + r - Thus, the arrival timing error is proportional to r, see fig. El 
(c). The steps correspond to the beating frequency which in general depends on f and 
f±. In this case when f = /_, it is (/+ — /_)/2 ~ 12 Hz. We note that this beating 
is crucial while fixing the coincidence timing window while performing coincidences 
between say Explorer- Nautilus or Explorer- Virgo. 

In case (b), when / is far away from FB, the situation is contrary. In this case, 
the timing error is small, [see Fig. 121(c)]. As explained earlier, for f away from FB and 
low r, the signal acts as a delta-like signal. As a result, delta-filter itself is a matched 
filter hence gives no timing error. Mathematically, the signal term in Eq. (jll|) is small 
compared to other terms. 
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Figure 3. Normalized output of the (A) matched filter, (B) delta filter, for t = 
IOtos, 0.2s, Is and /o = /_. 



A(a) tau = ID ms f =904.7 Hz A(bl tau = 200 ms f =904.7 Hz A(d tau = 1 s f =904.7 Hz 




6. Conclusion 

In this work, we did a comparative study of the response of the matched filter vs 
delta-filters for filtering the damped sinusoids GW signals for resonant bar detectors. 
We divided our study in two cases: (a) f G FB and (b) f far from the FB with 
FB = {900,932} using Explorer configuration. We find that in case (a), the loss in 
SNR increases as r and so does the arrival timing error if delta-filter is used instead of 
matched filter. However, in case (b), the signal almost acts like a delta for small r hence 
the SNR loss is negligible for small r however as r increases, the SNR loss gradually 
increases. The arrival timing error is minimal for all r. Thus, we can optimally use 
delta filters for detecting signals in case (b) for r as large as 200 ms. But, to detect 
damped sinusoids in case (a) with r as large as 50 ms, it is mandatory to use the optimal 
filtering as opposed to delta-filter. This poses a problem of setting an "optimal" grid 
of templates in f ,r when f ,r are fixed. However, as described earlier, f and r can 
evolve in the detector bandwidth. For detecting such signals and perform coincidences, 
the delta filter is inadequate and an optimal "matched" filter is difficult to construct due 
to insufficient knowledge of signal waveform. Alternative detection methods are needed 
which we pursue in future work. 
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